Phonon dispersion and electron-phonon interaction for YBa 2 Cu 3 07 from 

first-principles calculations 
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We present a first principles investigation of the lattice dynamics and electron-phonon coupling of 
the high-T c superconductor YBa2Cu307 within the framework of density functional perturbation 
theory using a mixed-basis pseudopotential method. The calculated phonon dispersion curves are in 
excellent agreement with Raman, infrared and neutron data. Calculation of the Eliashberg function 
a 2 F leads to a small electron-phonon coupling A=0.27 in disagreement with earlier approximate 
treatments. Our calculations strongly support the view that conventional electron-phonon coupling 
is not an important contribution to superconductivity in high-T c materials. 

PACS numbers: 74.24.Kc, 74.72.Bk, 78.30.-j 



Since the discovery of superconductivity in copper ox- 
ides in 1986 many different mechanisms have been sug- 
gested to explain the high transition temperatures ob- 
served in these materials jjj. However, even after more 
than 15 years of intensive research no consensus has been 
reached apart from the fact that conventional electron- 
phonon coupling alone cannot explain the observed tran- 
sition temperatures^, ^, ^] . Even the much simpler ques- 
tion concerning the magnitude of electron-phonon cou- 
pling has not been answered conclusively due to the fact 
that calculations in the past treated the lattice dynamics 
only in a very approximate way. Recent progress in the 
ab initio calculation of the lattice dynamics has made it 
possible to overcome these limitations. 

Modern bandstructure calculations together with ab- 
initio determination of the phonon dispersion, electron- 
phonon coupling and solution of the Eliashberg equations 
have made it possible to calculate transition tempera- 
tures and study in detail the electron-phonon coupling. 
This has recently been done very successfully for the 
newly discovered superconductor MgBg^. For high-T c 
superconductors (FTTSC) in the past however only a few 
attempts have been made to proceed along the above 
mentioned road map to a deeper understanding of the 
electron-phonon coupling M Et, ||, |?], ||, ||. The most 
involved ones are those concerned with La2Cu04 and 
doped CaCuC>2. For La2Cu04 the complete phonon dis- 
persion has been calculated however no electron-phonon 
coupling was determined. Complications arise from the 
fact that the calculations are carried out for the high-T c 
tetragonal phase while at low T c experimentally an or- 
thorhombic phase is observed. More important however 
is the fact that for undoped La2Cu04 an anitferromag- 
netic non-metallic ground state is observed while LDA 
calculations find a metallic ground state which questions 
the reliability of the calculations. The system CaCuC>2, 
which is much simpler, also requires hole doping which 
has been taken into account in an average way, however 
only after properly doping the Fermi surface becomes 
similar to those of other high-T c superconductors. Thus 
these results are also questionable. To our knowledge 



so far all other ab-initio treatments are restricted to the 
calculation of phonon modes for the P-point only. In 
addition, even these calculations are not complete but 
have been restricted to modes of special symmetry classes 
only. Due to these limitations estimates of the coupling 
strength are not very accurate, thus asking for a more 
complete treatment of the electron-phonon coupling in 
the whole Brillouin zone (BZ). 

To clarify the role of electron-phonon coupling in 
the copper oxides we have studied as a prototype 
YBa2Cu307 along the above mentioned lines. This sys- 
tem has been selected since it is the best experimen- 
tally studied high-T c superconductor with detailed in- 
formation available about the Fermi surface and the lat- 
tice dynamics and in addition does not have the com- 
plications of doping and/or structural phase transitions. 
Furthermore, bandstructure calculations based on the 
local-density approximation for exchange and correlation 
(LDA) have been successful in describing the Fermi sur- 
face and account properly for the metallic ground state 
0. For YBa2Cu307 we have carried out a systematic 
study of the lattice dynamics and the electron-phonon 
coupling using the mixed-basis pseudopotential method 
|fj"of . For all atoms we have used well-tested potentials of 
BHS-type [jn]. As valence states we included Y4s, Y4p, 
Ba5s, Ba5p and 02s which led to fairly deep potentials, 
but they could be dealt with very efficiently due to the 
mixed-basis formulation. The wave functions were con- 
structed from localized s,p,d functions at the yttrium and 
barium sites, s and p functions at the oxygen sites and 
d-functions at the copper sites supplemented by plane 
waves with an energy cut-off of 20 Ry. Detailed tests 
were carried out to ensure convergence with respect to 
the number of plane waves as well as with respect to 
k-point sampling. For the structural optimization 576 k- 
points in the BZ together with a Gaussian smearing of 
0.2 eV have been used. The calculation of the phonon 
dispersion is based on a recently developed mixed ba- 
sis perturbation approach jl^] which also allowed for the 
calculation of the electron-phonon coupling 13 . Since 



these parts are very k-point sensitive up to 5184 k-points 
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FIG. 1: Crystal structure of YBa 2 Cu 3 7 



TABLE I: Comparison of calculated and measured structural 
parameters of YBa2Cu3C>7. 



this work 



Vol (A a ) 
c/a 
b/a 

Zflo 
Zcu(2) 
Z 0(2) 
Z 0(3) 
Z 0(4) 



FLAPW[6| 



Exp. 15, 16] 



165.6 
3.015 
1.021 
0.1819 
0.3495 
0.3770 
0.3761 
0.1627 



163.1 
3.015 
1.012 
0.1817 
0.3507 
0.3770 
0.3765 
0.1619 



172.32 
3.056 
1.018 
0.1837 
0.3546 
0.3783 
0.3779 
0.1592 



in the BZ have been used here. All our calculations were 
carried out using the local exchange-correlation potential 
of Hedin and Lundqvist jl4j . 

Since YBa 2 Cu 3 07 crystallizes in the orthorhombic 
crystal structure optimization requires the determination 
of 8 parameters (the length of the axis a, b, c and 5 in- 
ternal positions) (see Fig. 1). In tabic I our results are 
summarized and compared to all- electron results || as 
well as experimental data 111], 16 . 



The two calculations give very similar results and the 
differences are well within the limit of numerical accu- 
racy. The calculated equilibrium volume turns out to be 
4% to 5% smaller than the experimental one in agree- 
ment with previous calculations ^] and general trends 
observed for LDA-calculations. 

Having determined the equilibrium structure we now 
proceed to the calculation of phonon modes. For compar- 
ison with other calculations as well as experimental stud- 
ies we concentrate first on the T point. Here, besides neu- 
tron scattering studies also Raman and infrared data are 
available. Table II summarizes the even symmetry modes 
and compares our results with those of an all-electron 
study ||(full potential linearized augmented plane- wave 
calculations (FLAPW)) which is to our knowledge the 
most extensive phonon study available up to now as well 
as with experimental data |l7], As pointed out in 



TABLE II: Calculated and measured Raman modes. Besides 
the frequency (in meV) also the dominant eigenvector char- 
acter is given. 



Mode 



this work 



FLAPW§ Exp.]l7j 



Ai 


Ba 


14.8 


15.3 


14.4-14.8 




Cu(2) 


18.8 


18.2 


18.0-18.6 




0(2)-0(3) 


42.1 


41.9 


41.5 




0(2)+0(3) 


50.4 


52.3 


53.9-54.6 




0(4) 


58.6 


60.4 


61.1-62.0 


B 2s 


Ba 


8.0 


8.1 


8.7 




Cu(2) 


17.4 


17.6 


17.6 




0(4) 


27.8 


27.5 


26.0 




0(3) 


47.9 


48.2 


45.9 




0(2) 


71.4 


73.5 


71.8 


B 39 


Ba 


9.8 


9.8 


10.3 




Cu(2) 


17.4 


17.5 


17.4 




0(4) 


36.6 


36.3 


37.6 




0(2) 


45.3 


46.1 


46.9 




0(3) 


65.2 


67.7 


65.2 



TABLE III: Calculated and measured odd symmetry modes 
for the r point in meV. 





Mode 


Theory 


Neutron 


IR 


Bi u 


Ba 


16.1 








CJu(l) 


20.7 


19.4 


19.2 




Y/Cu(2) 


23.8 


23.9 


23.9 




0(2)-0(3) 


27.0 




silent 




0(1) 


29.8 




34.1 




0(2)+0(3) 


40.3 


37.2 


37.2 




0(4) 


68.2 


68.1 


70.3 


B2u 


Y/Ba 


10.3 


10.2 






0(1) 


16.4 


15.6 






Cu(l) 


19.0 


19.5 






Y/Cu(2) 


21.8 


24.0 


23.5 




0(4) 


33.5 


35.0 


34.7 




0(3) 


44.6 


44.0 


44.9 




0(2) 


71.9 


74.3 


74.0 




Y/Ba 


10.0 


10.1 






Cu(l) 


19.7 


19.5 






Y/Cu(2) 


21.2 


23.4 


23.4 




0(4) 


41.7 




42.7 




0(2) 


42.6 


44.5 


44.9 




0(1) 


61.4 


59.4 


59.2 




0(3) 


65.6 


68.0 


67.6 



connection with the FLAPW-study structure optimiza- 
tion has an important effect on phonon frequencies. This 
has also been found in our study. The agreement between 
both theoretical treatments and the Raman data is very 
satisfactory. The pseudopotential treatment is of equal 
quality as the FLAPW-method, in addition it is compu- 
tationally less demanding and allows thus the calculation 
of odd modes as well which are given in table III together 
with neutron- and IR-data |l^, |T^] . 

So far, we have discussed T-point modes only, how- 
ever, a reliable calculation of electron-phonon coupling 
requires the knowledge of the complete phonon disper- 
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FIG. 2: Phonon dispersion along [100] direction. The modes 
have been classified according to symmetry. Lines represent 
theoretical results while dots are neutron data 1181. 
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FIG. 3: Phonon dispersion along [010] direction. The modes 
have been classified according to symmetry. Lines represent 
theoretical results while dots are neutron data Il8|. 



sion throughout the BZ. Results for different directions in 
the BZ obtained using our mixed basis perturbation ap- 
proach (l2| are presented in Figs. 2 and 3. The modes are 
classified according to symmetry. For comparison also 
experimental data points are indicated. The dispersion 
curves were obtained by determining the dynamical ma- 
trix on a (4,4,2) reciprocal lattice grid in the BZ. Keeping 
in mind the complexity of the system the results are in 
excellent agreement with the available neutron data. 

A number of points are worth mentioning in connection 
with the dispersion curves. The strong renormalization 
in respect to YBa2Cu3C>6 of the top- most mode of Ai- 
symmetry along [100] shows up quite naturally in these 
calculations while model studies had extreme difficulties 
in reproducing this behavior. The eigenvector shows that 
it is an oxygen bondstretching mode in the Cu-0 plane. 




20 40 60 

CO (meV) 



FIG. 4: Upper panel: Generalized phonon density of states. 
A broadening of 0.5 meV has been used. Lower panel: Eliash- 
berg function a 2 F(ui) (full line) and phonon density of states 
F(w) (dotted line). 



For experimentalists one of the major obstacles in the 
past has been the lack of untwinned single crystals of 
YBa2Cu3C>7 which made sometimes the assignment of 
a measured phonon to the (100) or (OlO)-direction very 
difficult. Some of these ambiguous cases have been re- 
assigned in view of the calculations and in agreement 
with very recent experimental studies [p0| . Concerning 
the role of the Cu-0 chains we found that most of these 
modes are in the range between 10 and 20 meV and only 
oxygen modes along the bond direction are high-lying 
ones with frequencies around 60 meV. No indication of a 
possible instability of the bond-bending chain modes (as 
proposed by earlier calculations [^))was found in agree- 
ment with experiment. Consistent with the expectation 
of weak coupling perpendicular to the Cu-O-planes the 
dispersion in (001) direction is very small. 

In the upper panel of Fig. 4 we show the generalized 
phonon density-of-states which should be compared with 
the neutron-scattering results presented in Refs. ^l], p2[ . 
Peaks and even shoulder positions agree quite well, as 
also the frequency cut off of 75 to 80 meV. The ratio 
of peak intensities differs substantially between the two 
measurements which is most likely due to the sensitivity 
of intensities to the background corrections. This pro- 
hibits comparison with our calculation. 

To address the superconducting properties we have cal- 
culated the so-called Eliashberg function a 2 F(uj) us- 
ing a very fine (24,24,4) k-point grid in performing the 
Fermi-surface average of the electron-phonon matrix el- 
ements j24j. Extensive tests have shown that this grid 
is sufficient for our calculations. Results are presented 
in Fig. 4 (lower panel). Comparing the phonon density 
of states F(oj) and a 2 F(ui) we find a substantial shift of 
weight for a 2 F(w) to higher frequencies. 
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From the Eliashberg function we calculate the electron 
phonon coupling constant A = 2 L duja 2 F (uj) / uj which 
gives A=0.27. For the logarithmically averaged phonon 
frequencies as defined in Ref. |E5| we find 27.7 meV. Af- 
ter solving the isotropic gap equation for ^*=0 we obtain 
only a T c of ~ 2 K. Any realistic /i* suppresses super- 
conductivity completely. 

The very small value of A is in substantial disagreement 
with earlier estimates based on T-point results only. An- 
alyzing the behavior of A throughout the BZ we find it 
to be varying by a factor of two only, however, varia- 
tion as function of individual modes for given q is much 
larger. For the T-point for example modes of A g and 
Bi u -symmetry account for 80% of A. The largest con- 
tributions come from modes in the copper-oxygen planes 
as expected. Previous estimates of A assumed that all 
modes at T contribute with similar strength as the A g - 
modes. This elucidates why A-values of the order of 1 or 
even larger have been obtained. However, this is a pure 
artefact of the unjustified assumption. 

Raman measurements at T allowed for the determi- 
nation of the anomalous frequency shifts of lattice vi- 
brations occurring when cooling a superconductor below 
T c p6[ . Based on a theory by Zeyher and Zwicknagl 
for the phonon-self energy it could be shown that val- 
ues of Ai=0.02 for the A g mode around 42 meV and 
A2=0.01 for the A ff -mode around 54 meV, as determined 
by an LMTO-calculation @, are consistent with the ob- 
served experimental frequency shifts. Our calculation 
give slightly larger values of Ai «0.06 and A2 «0.02 which 
are equally consistent with the measurements due to the 



approximations inherent in the calculation of the phonon- 
self energy [^7| . 

So far we have discussed only s-wave pairing, however 
many experimental facts support d-wave pairing |]28| . In 
the much simpler system CaCu02 the influence of s- ver- 
sus d-pairing has been studied [Q and it has been found 
that both channels give very similar results for A which 
should also hold for YBa2Cu30y. One other question 
which has not been discussed so far is the influence of 
correlations beyond LDA. Here the strongest argument 
in favor of our calculations is the fact that the calculated 
phonon dispersion is in very good agreement with exper- 
iment. If the clcctron-phonon coupling would be very 
strongly dependent on correlations this should show up 
in significant discrepancies between calculated and mea- 
sured phonon modes. The good agreement supports the 
claim that correlation effects are of minor importance for 
electron-phonon coupling. 

In summary, we have presented here ab-initio results 
for the lattice dynamics of YBa2Cu30y which are in ex- 
cellent agreement with available Raman, infrared and 
neutron data. Based on the knowledge of the complete 
phonon spectrum the electron-phonon coupling has been 
calculated. The electron-phonon coupling constant turns 
out to be very small. These calculations strongly support 
the view that in high-T c materials conventional electron- 
phonon coupling is not an important contribution to su- 
perconductivity. 

The authors thank Dr. L. Pintschovius and Dr. W. 
Reichardt for numerous discussions. 
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